%*******************************************************************************
%*
%*   Copyright ©2009 Alexandru Marin
%*
%*   This file is part of nwave/pMRI C library.
%*
%*   This project started as an MSc degree project during a
%*   stage at the Paris-Est university, under joint supervision of prof.
%*   Bogdan Ionescu (Image Analysis and Processing Laboratory, "Politehica"
%*   University of Bucharest), Caroline Chaux and Jean-Christophe Pesquet
%*   (Laboratoire d'Informatique, Université Paris-Est, Institut Gaspard Monge).
%*
%*   nwave/pMRI C libray is free software: you can redistribute it and/or modify
%*   it under the terms of the GNU General Public License as published by
%*   the Free Software Foundation, either version 2 of the License, or
%*   (at your option) any later version.
%*
%*   nwave/pMRI C libray is distributed in the hope that it will be useful,
%*   but WITHOUT ANY WARRANTY; without even the implied warranty of
%*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%*   GNU General Public License for more details.
%*
%*   You should have received a copy of the GNU General Public License
%*   along with nwave/pMRI C library. If not, see <http://www.gnu.org/licenses/>
%*
%*******************************************************************************

%
%
% This example ilustrates the results of the DWT/IDWT 
% from the mrilib package
%
% You first need to compile and run the dwt_example.c file
%

original = mri_read_matrix('testimage.csv');
figure;image(original);colormap(gray(256));

transformed = mri_read_matrix('testimage_transf.csv');
figure;image(transformed);

reconstructed = mri_read_matrix('testimage_rec.csv');
figure;image(reconstructed);colormap(gray(256));

disp(['psnr versus original image: '...
    num2str(mri_psnr(original, reconstructed))]);